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Abstract. We propose a simple schematic model for two-neutron halo nuclei. In 
this model, the two valence neutrons move in a one-dimensional mean field, interacting 
with each other via a density-dependent contact interaction. We first investigate the 
ground state properties, and demonstrate that the dineutron correlation can be realized 
with this simple model due to the admixture of even- and odd-parity single-particle 
states. We then solve the time-dependent two-particle Schrodinger equation under the 
influence of a time-dependent one-body external field, in order to discuss the effect of 
dineutron correlation on nuclear breakup processes. The time evolution of two-particle 
density shows that the dineutron correlation enhances the total breakup probability, 
especially for the two-neutron breakup process, in which both the valence neutrons 
are promoted to continuum scattering states. We find that the interaction between 
the two particles definitely favours a spatial correlation of the two outgoing particles, 
which are mainly emitted in the same direction. 
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1. Introduction 

Neutron- rich nuclei have attracted much interest during the past decades [U Ej |3l HI El E] , 
and this will continue to be so due to new generation radioactive beam facilities in 
the world. These nuclei are characterized by a small binding energy, and many new 
features originating from the weakly bound property have been found. A halo and 
skin structures with a large spatial extension of the density distribution [7], a narrow 
momentum distribution [8j, a modification of shell structure and magic numbers^, 
and strong concentration of electric dipole {El) transition at low excitation energies 
[TOl im [12] . are well-known examples. 
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Among neutron-rich nuclei, two-neutron halo nuclei are particularly intriguing 
systems to study. Their structure has often been described as a three-body system 
consisting of two valence neutrons interacting with each other and with the core nucleus 
[U |5l [131 [H]. Some light neutron-rich nuclei, such as ^^Li and ^He, do not have a bound 
state in the two-body subsystem with a valence neutron and a core nucleus. These nuclei 
are referred to as Borromean nuclei, and their properties have been studied extensively 
both experimentally [B 13 [3 El HI |I2] and theoretically [3lSl|5l[6l[l5l|16l|l71ll8]. Nuclear 
breakup reactions of Borromean nuclei have also been investigated by the continuum- 
discretized-coupled-channels (CDCC) method [HI [201 [211 [221 [23] and by the eikonal 
method |2il [251 [26] . 

One of the most important current open questions concerning the Borromean nuclei 
is to clarify the characteristic nature of correlations between the two valence neutrons, 
which do not form a bound state in the vacuum. A strong dineutron correlation, 
where the two neutrons take a spatially compact configuration, has been theoretically 
predicted[3l [H [5l [161 EH [281 [SSI ES]- Although the recent experimental observation of 
the strong low-lying dipole strength distribution in the ^^Li nucleus [12] has provided 
an experimental signature of the existence of dineutron correlation in this nucleus, its 
direct evidence has not yet been obtained. 

The pair transfer reaction is another promising way to probe the dineutron 
correlation, as the cross section is known to be sensitive to the pairing correlation [3T1 [32] . 
However, the reaction dynamics is rather complicated and the role of dineutron in the 
pair transfer reaction has not yet been fully clarified, although experimental studies 
on pair transfer reaction with exotic nuclei have been initiated recently [331 [311 [35] . 
Apparently, it is urged to construct a theoretical framework for pair transfer which 
fully takes into account the pairing correlation and the dineutron correlation in its 
consequence. 

The aim of this paper is to develop a simple toy model for two-neutron halo nuclei. 
Albeit its schematic nature, such model is rather useful as it allows detailed studies 
on the static and dynamical properties of two-neutron halo nuclei. A schematic model 
is also useful to deepen our understanding of two-neutron halo nuclei by providing 
intuitive pictures of several dynamical processes. To this end, we consider a three-body 
model in which the motion of two valence neutrons is confined within one-dimensional 
spatial space. The neutrons are assumed to move in a one- dimensional Woods-Saxon 
potential, while interacting with each other via a two-body interaction, which we take 
a density-dependent contact interaction [H [151 US]. This is a natural extension of the 
model developed in Ref. [36] for a one-neutron halo nucleus. We apply the model to the 
ground state as well as to a nuclear breakup process of a loosely-bound two-neutron halo 
nucleus, leaving an application to the pair transfer process in a separate publication. In 
passing, a similar one-dimensional three-body model has been used in atomic physics, 
in order to discuss e.g., the mechanism of the double ionization process of He atom by 
intense laser fields [371 EH [39] . 

The paper is organized as follows. In Sec. II, we detail the one- dimensional three- 
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body model. In Sec. Ill, we apply the model to the ground state of a loosely-bound 
two-neutron halo nucleus and discuss the dineutron correlation. In Sec. IV, we consider 
the two-neutron halo nucleus under the influence of a time-dependent external fleld. We 
solve the time-dependent Schrodinger equation in one-dimensional coordinate space and 
discuss the time evolution of the two-particle wave function. We then summarize the 
paper in Sec. V. 



2. One-dimensional three-body model 

We consider a two-neutron halo nucleus in a one dimensional space. Denoting the 
coordinate of the two neutrons by xi and X2, the three-body Hamiltonian reads, 

H = + V{x,) - + V{x,) + vUx., X,), (1) 

where m is the nucleon mass. Here, we neglect for simplicity the recoil kinetic energy of 
the core nucleus. V{x) is the neutron-core potential, which we take the Woods-Saxon 
form [36], 



Vnn is the neutron-neutron interaction. We take the density- dependent contact 
interaction for it, that is, 

Vnn{x, X') = -g(l- Yl^^^mr^ ^("^ ~ 

where we assume that the density is given by the same functional form as the mean-fleld 
potential, Eq. ([2]). This is the so called surface type pairing interaction, which almost 
vanishes near the center of the core nucleus at x ~ 0. The density- dependent contact 
interaction has been successfully employed in the Hartree-Fock-Bogoliubov calculations 
for neutron- rich medium-heavy and heavy nuclei |28l HO] . as well as in describing the 
structure of ^^Li and ^He nuclei with the three-body model [H [151 116] . 

In order to obtain the ground state wave function of the two-neutron halo nucleus, 
we flrst solve the Schrodinger equation for the two-body subsystem, 

2m dx"^ 

Each eigenstate n is assumed to have a two-fold degeneracy, depending on the direction 
of the spin of neutron. The continuum states are discretized by putting the nucleus 
within a large box. Denoting the size of the box to be Xbox, we impose the vanishing 
boundary condition 0„(x) = at x = ±Xbox- 

We expand the two-particle wave function \E'(xi,X2) with the single-particle wave 
functions (j)n{x) as 



+ V{x) - e„ 



<Pn{x) = 0. (4) 



^gs(xi,X2) = ^ a„„/^„„/(xi,X2), (5) 

n<n' 



Two-neutron halo nuclei in one dimension 



4 



where 

'^nn'{Xi,X2) = , ^ = [(j)n{Xl)(j)n'{x2) + 4>n{x2)(j)n'{Xl)] \S = 0). (6) 
.^2(1 + Sn,n') 

Here we have assumed that the spin of the two neutrons form the total spin of S=0 
configuration, so that the wave function is symmetric with respect to the interchange 
of Xi and X2- Notice that the Hamiltonian ([T]) conserves the parity. Since the ground 
state has positive parity, the expansion in Eq. ([5]) can therefore be restricted to the 
configurations where the states n and n' have the same parity. The expansion coefficients 
ann' are determined by diagonahzing the Hamiltonian matrix, whose matrix elements 
read. 



2(l + 5ni,nJv/2(l + 5„;xJ 

(7) 

Once the ground state wave function \E'gs(a;i, X2) is obtained, the two-particle and 
one-particle densities can be constructed as 

P2(a;i,X2) = |^gs(a;i,X2)P, (8) 

and 

Pi{x) = J dx'p2{x,x'), (9) 

respectively. The mean value of the neutron-neutron distance and the distance between 
the core and the center of mass of the two neutrons are also computed as 



Xnn =\ dXidX2{Xi - X2yp2{Xi,X2), (10) 



fXi + X2^ ^ 



Xc2n = \ll dXidX2 ) p2{Xl,X2), (U) 



respectively. 



3. Ground state properties 

We now apply the one- dimensional model to a weakly-bound two-neutron halo nucleus. 
To this end, we take R = 1.27 x 9^/^ fj^^ ^ ^ -50.085 MeV, and a=0.67 fm for the 
Woods-Saxon potential, Eq. ([2]), in order to mimic the ^^Li nucleus. This potential 
possesses four bound single-particle states. We assume that the lowest three bound 
states are already occupied by the neutrons in the core nucleus, and exclude those 
states explicitly from the expansion in Eq. (jS]). In this way, the Woods-Saxon potential 
has effectively only one bound state with odd parity at ews = —0.15 MeV. That is, the 
three-body system is bound by 0.3 MeV without the pairing interaction. Although we 
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could set up a situation in which there is no bound state in a Woods-Saxon potential, 
we introduce a loosely bound single-particle state in order to mock up a resonance state, 
which plays an important role in the ^^Li and ^He nuclei [T7] . 
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Figure 1. The dependence of the strength of the pairing interaction g on the cutoff 
energy, i?cut- For each i?cut, the strength g is adjusted in order to yield the ground 
state at E^.s. = -1 MeV. 



It is well known that a contact interaction has to be supplemented with an energy 
cutoff, Ecut in US]- For a three-dimensional case, the strength of the pairing interaction, 
g, is related to the cutoff energy, -Ecut, via a scattering length a^n for the n + n scattering 
[ll[15]. On the other hand, for a one- dimensional case, we simply vary the value of g for 
each cutoff energy (with keeping the Woods-Saxon potential well) in order to reproduce 
a given ground state energy, since it is not straightforward to define the phase shift for 
a one-dimensional scattering problem. Figure 1 shows the dependence of g on Ecut for 
the ground state energy of E^s. = — 1 MeV, obtained by taking the configurations in 
Eq. ([5]) which satisfy + e„/ < Ecut- One notices that the dependence is rather weak. 
In this paper, we arbitrarily take the cutoff energy to be 30 MeV, that corresponds to 



-12.35 MeV- fm for E, 



g.s. 
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Figure 2. The energy spectrum for the model Hamiltonian used in this paper. The 
energies are measured from the threshold of a three-body scattering state. The three- 
body bound state, [c+n+n], is located at i? = — 1 MeV. The [n + n] and [c + n] 
are two-body bound states for the two-neutron (n-l-n) and the core+neutron (c-|-n) 
systems, respectively. 



It is important to notice that a one-dimensional delta function potential v{x) = 



Two-neutron halo nuclei in one dimension 



6 



—g6{x) always holds a bound state at E = —mg'^/Ah'^ for a two-neutron system even 
with an infinitesimally small attraction g For our choice of g = —12.35 MeV- fm, 
a dineutron is thus bound by 0.92 MeV. See Fig. 2 for a spectrum for the three-body 
system considered in this paper. There is no bound excited state below the threshold of 
the core nucleus -|- a bound dineutron ai E = —0.92 MeV. Thus, the ground state in the 
present model has an extremely small separation energy from the threshold for the di- 
neutron breakup, i.e., less than 100 keV. Since the dineutron is bound, our calculations 
correspond to a limit of strong neutron-neutron interaction in a three-dimensional case. 
Alternatively, our calculations also have a similarity to ^Li in which a bound deuteron 
may exist in the channel of ^Li = a + p + n. 
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Figure 3. The energy of the last bound single-particle level (Fig. 3(a)), the ground 
state energy of the three-body system (Fig. 3(b)), the neutron-neutron mean distance 
(Fig. 3(c)), and the core-dineutron mean distance (Fig. 3(d)) as a function of the size 
of the box, A'box- 



Figure 3 shows several quantities for the ground state of the three-body 
Hamiltonian, Eq. ([1]), as a function of the size of the box, Xbox- Figs. 3(a), 3(b), 
3(c), and 3(d) are for the single-particle energy of the last bound level, the ground 
state energy of the three-body system, the neutron-neutron root mean square distance, 
and the mean distance between the core and the center of mass of the two neutrons, 
respectively. These quantities are almost converged at Xbox ~ 50 fm. For a Borromean 
system, it was shown that the convergence of the ground state wave function is much 
slower than the ground state energy, especially in the tail region [l2lll3]. For the example 
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shown in Fig. 3, in which there is a bound single-particle state, we have confirmed that 
the wave function has also converged at Xbox ~ 50 fm. In the calculations shown in this 
paper, we take the box size to be Xbox=90 fm. 
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Figure 4. (Color online) The two-particle density for the correlated (Fig. 4(a)) and 
the non-correlated (Fig. 4(b)) ground states of a one-dimensional three-body model. 



Figure 4(a) shows the two-particle density of the ground state. The density 
distribution is largely concentrated along the line of xi ^ X2, that is nothing but the 
manifestation of strong dineutron correlation discussed in Refs. |ll O [161 1211 [28l [291 130] ■ 
The correlation largely hinders the density along the Xi = —X2 line, and only the peaks 
along xi ~ X2 survive. When the pairing interaction is switched off, the two-particle 
density becomes that shown in Fig. 4(b). The density for the non-correlated ground 
state is symmetric with respect to the transformation of xi — > —xi {x2 — )■ —X2) for a 
fixed value of X2 {xi), and it has therefore four symmetric peaks. See Ref. [l^ for a 
similar figure for the ^^Li nucleus. 



4.0x10"'* (b) 20 



3.0x10" 



10 



2.0x1 0"'* ^ 
1.0x10"'* -10 



0.0x10 



,0 



-20 





•1 




•1 



0.007 (c) 20 
0.006 

0.005 ^ 10 
0.004 1, n U 
0.003 CM 
0.002 ^ .-[Q 
0.001 

-20 



-20 -10 10 20 
Xi (fm) 




0.002 
0.001 


-0.001 
-0.002 



-20 -10 10 20 
xi (fm) 



-20 -10 10 20 
xi (fm) 



Figure 5. (Color online) The decomposition of the correlated two-particle density 
into the even-parity contribution (Fig. 5(a)), the odd-parity contribution (Fig. 5(b)), 
and the interference between the even and odd parity contributions (Fig. 5(c)). 



It is well known that the dineutron correlation is caused by the admixture of many 
single-particle states with different parity [291 [15]. order to demonstrate it explicitly, 
we decompose the ground state wave function, Eq. ([5]), into two components, 

^fgs(xi,X2) = ^ee(a;i,a;2) + ^'oo(a;i,a;2), (12) 
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where \l/ee(a;i, 2:2) and \l/oo(a^i, 2^2) consist only of the even-parity and the odd-parity 
single-particle states, respectively. Assuming that \E'ee and \&oo are both real, the two- 
particle density then reads. 

Fig. 5 shows each of the components separately (see also Fig. 12 in Ref. [29j). With 
only even or odd parity single-particle states, the density distribution has four symmetric 
peaks as in the non-correlated case shown in Fig. 4(b). l^&ooP is much more compact 
than I^I/eeP because the former contains the bound single-particle state. It is clearly seen 
in Fig. 5(c) that the localization of the two- neutron pair along the line Xi = X2 emerges 
only with the interference between the even and odd parity states [29ll45]. 

This can be understood easily as follows. Because the nuclear interaction is short 
ranged, the main contribution comes from n = n' in Eq. ([5]). Suppose that there are 
only one even-parity single-particle state 4>e{x) and one odd-parity state, 4>o{x). If we 
write the ground state wave function as 

\E'g.s.(a;, x') = a (f)e{x)4>e{x') + P (j)o{x)(j)o{x'), (14) 
the two-particle density reads 

P2{x,x') = |0e(a;)0e(a;')P + /^^ \4>o{x)4>o{x')'^ 

+ 2a/3 0e(a;)0e(x')0o(a;)0o(a;')- (15) 

As we prove in the Appendix by using the two-level model, a and /3 have the same sign 
for the density-dependent pairing interaction. The first and the second terms in this 
equation are positive definite. On the other hand, the third term is positive for x = x', 
while it is negative for x = —x' as (j)e{—x) = (f)e{x) and (j)o{—x) = —(f)o{x). Therefore, the 
interference term is destructive and decreases the two-body density along the x = —x' 
line, while it is constructive and enhances the density along the x = x' line. 

In the past, the dineutron correlation was discussed using several representations. 
In Refs. [H |28l Hi] , the density distribution of the second neutron was investigated when 
the first neutron was fixed at a certain position. In Refs. [5|[29H27t[l5t[36l[i71[i8Hl9], 
the two-particle density was plotted as a function of the relative distance between the 
neutrons, r = ri — r2, and the center of mass coordinate, R = [vi + r2)/2. In Ref. 
[T6] . by setting ri = r2 = r the two-particle density was plotted as a function of 
the core- neutron distance, r, and the opening angle between the two neutrons, 6*12. 
Although all the representations are useful to reveal the strong dineutron correlation, 
they have advantages and disadvantages. For instance, it is not easy to explore all 
the position of the first neutron in the first representation, and the configurations with 
ri 7^ r2 are neglected in the third representation. These can be avoided with the second 
representation, but it is more intuitive to use directly the coordinates ri and r2, rather 
than r and R, especially when one has to consider also the angular part of r and R. The 
one-dimensional model removes these inconveniences, yielding a transparent picture for 
the two-particle density distribution. 
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4. Nuclear Breakup process 

Let us next discuss the time evolution of the two-particle wave function in the presence 
of an external field acting on each particle. As the external field, we take [36] 

Kxt(xi,X2,t) = E v;e-*'/2'^?e-(^-^»)'/2.^, (16) 

2 = 1,2 

with the parameters of V^=3 MeV, at = 2.1h/MeV, and ax=2.2 fm. In order to 
investigate the time evolution, we solve the time-dependent two-particle Schrodinger 
equation, 

d 

ih—^{xi,X2, t) = [H + V;xt(a;i, X2, t)]^(xi, xs, t), (17) 
with the initial condition of 

^(a;i, X2, to) = *gs(a;i, 3:2), (18) 

at the initial time to. We ortho normalize the wave function to the occupied states 
using the projection procedure at each time step during the time evolution. In the 
calculations shown below, we take cto = —400 fm, and use the implicit method [50j for 
the time propagation. We consider the symmetric perturbation, that is, = in Eq. 
f fT6|l . We have confirmed that our conclusions remain qualitatively the same (except for 
the asymmetry in the two-particle density distribution along the xi = X2 line) even if 
we choose an asymmetric perturbation [36j, e.g., with xq = 2 fm. 

Figures 6(a), 6(b), 6(c), and 6(d) show the two-particle density at ct = —400, 0, 300, 
and 600 fm, respectively. As the time evolves, the extension of the peaks along the 
Xi = X2 line increases significantly. This is in marked contrast to the uncorrelated 
case shown in Figs. 7(a), 7(b), 7(c), and 7(d). In the uncorrelated case, the two-particle 
density expands democratically, indicating that there is the equal probabihty of emission 
of the two neutrons in the opposite directions to that in the same direction. 

The corresponding current distributions, 

ji{Xi,X2,t) = -— {'$*{xi,X2,t)—'${xi,X2,t) 
2X171 \ OXi 

- ^{XuX2,t)—<i/*{x,,X2,t)\ , (19) 

are shown in Figs. 6(e-g) and 7(e-g) for the correlated and the uncorrelated cases, 
respectively. For the correlated case, the main flows are the outgoing flows along the 
Xi = X2 line, as one can infer from Figs. 6(b), 6(c), and 6(d). For the uncorrelated 
case, on the other hand, there are four symmetric outgoing flows that correspond to the 
expanding density distribution shown in Figs. 7(b), 7(c), and 7(d). 

In order to see more clearly the time evolution of the breakup fragments, we project 
the two-particle wave function at each time t onto the component which is orthogonal 
to the ground state. That is, 

^bu(a;i, X2, t) = ^(xi, X2, t) - (*gs|^(t)) ^gs(xi, X2). (20) 
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Figure 6. (Color online) The time evolution of the two-particle density. Figs. 6(a), 
6(b), 6(c), and 6(d) show the two-particle density, |^'(a;i, a;2)P at ct = —400,0,300, 
and 600 fm, respectively. Figs. 6(e), 6(f), and 6(g) show the corresponding current 
distributions with arbitrary units. Figs. 6(h), 6(i), and 6(j) show the breakup 
component of the two-particle density, |5'bu(a;i, a;2)P at ct — 0,300, and 600 fm, 
respectively. Notice the difference scales among Figs. 6(a), 6(b), 6(c), and 6(d). 
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Figure 7. (Color online) Same as Fig. 6, but for the uncorrelated case with a vanishing 
neutron- neutron interaction. Notice the difference scales among Figs. 7(a), 7(b), 7(c), 
and 7(d). 
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Figure 8. (Color online) The decomposition of the breakup component of two-particle 
density at ct = 600 fm shown in Fig. 6(j). Fig. 8(a) shows the (be) component, in 
which one of the neutrons is in a continuum state while the other remains in the 
bound single-particle state. Fig. 8(b) shows the (cc) component, in which both the 
two neutrons are in continuum states. 



Since there is only one three-body bound state (the ground state) below the dineutron 
threshold, all the remaining wave functions fl2U]) necessarily correspond to the breakup 
processes. Figs. 6(h), 6(i), and 6(j) show the time evolution of the breakup component 
of the two-particle density, |\l/bu(a;i, X2,t)\^. One can clearly see that the two neutrons fly 
apart from the core nucleus by sticking to each other due to the final state interaction. 
The effect of the correlation is clearly evidenced by the comparison with the uncorrelated 
case, shown in Figs. 7(h), 7(i), and 7(j). 

In order to clarify the dynamics of breakup process, we further decompose the 
breakup component of the two-particle wave function into three contributions, 

^bu(a:i, X2, t) = X2, t) + ^t"\xu X2, t) + ^(^'^^(Xi, X2, t), (21) 

where \l/bb and \E'cc correspond to the component in which both the neutrons are in the 
bound single-particle state and in the continuum states, respectively. \l/bc describes the 
component in which one of the neutrons is in the bound state while the other is in a 
continuum state. That is, 

^t^(xi,X2,t) = ^00(X1,X2), (22) 

^t)(xi,X2,t) = E P^n''\t) ^on(a;i,X2), (23) 

Vl/(^)(Xi,X2,t)= E 72^(t)^nn'(^l,^2), (24) 

n<n' 
{n,nVO) 

where '^nn' is given by Eq. (j6]), n = corresponding to the bound single-particle level 
at —0.15 MeV. Figs. 8(a) and 8(b) show the (be) and (cc) components of the breakup 
density, |^'bc(3^i, a;2)P and |\l/cc(a;i, X2)P, respectively, at ct = 600 fm shown in Fig. 6(j). 
The (bb) component is small, and is not shown in the figure (the norm of \E'bb, ^bc, 
and \l/cc are 4.58x10"'*, 0.116, and 0.466, respectively. See also Fig. 11 below). It is 
interesting to notice that the inner part of the density for the (be) component somewhat 
resembles the non-correlated density at early times shown in Fig. 7(h). As we will show 
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in Fig. 10, this originates from the fact that the (be) component makes a dominant 
contribution in the uncorrelated case. For the (cc) component, on the other hand, the 
two emitted neutrons are sticking to each other due to the final state interaction, and 
the breakup process may be interpreted as an emission of a bound dineutron. 
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Figure 9. (Color online) The two-particle densities at ct = 600 fni obtained by 
neglecting the final state interaction (i.e., the neutron- neutron interaction) during 
the time evolution from the correlated ground state at initial time. Figs. 9(a) and 
9(b) show the total density, \'i'{xi,X2)\'^ , and its breakup component, |^'bu(a;i, a;2)|^, 
respectively, while Figs. 9(c) and 9(d) show the (be) and the (cc) components of the 
breakup density, respectively. 



We point out that both the dineutron correlation in the ground state, as well as 
the neutron-neutron interaction acting during the time evolution (that is, the final state 
interaction) are important for the dineutron emission process. In order to demonstrate 
this. Fig. 9 shows the two-particle density at ct = 600 fm obtained by switching off 
the neutron-neutron interaction during the time evolution. The initial state is still 
prepared as the correlated ground state. Figs. 9(a) and 9(b) show the total density, 
|\E'(a;i,a;2)P, and its breakup component, |^'bu(a;i, a;2)P, respectively. The (be) and 
the (cc) components of the breakup density are also shown in Figs. 9(c) and 9(d), 
respectively. One can see that the density distributions are considerably different from 
those corresponding to the fully correlated case. Figs. 6(d), 6(j), 8(a), and 8(b). Without 
the final state interaction, the probability for the non-correlated two-neutron emission 
is increased significantly (see the peaks in Figs. 9(a), 9(b) and 9(d) along the Xi = —X2 
line). At the same time, the dineutron emission is somewhat slowed down, as the peaks 
in Figs. 9(d) along the Xi = X2 line are located much closer to the core nucleus as 
compared to the fully correlated case shown in Fig. 8(b). 
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Figure 10. (Color online) Same as Figs. 9(a) and 9(b), but for the case with the 
final state interaction during the time evolution from the uncorrelated ground state at 
initial time. 



Figure 10(a) and 10(b) show the calculated results starting with the uncorrelated 
initial state, but with taking into account the neutron-neutron interaction during the 
time evolution. To this end, we multiply the same time profile function, e"*^/^'^' , as in the 
time-depenent external field, Eq. [161 to the pairing interaction in the time-dependent 
Schrodinger equation. Figs. 10(a) and 10(b) show the total density, |\E'(xi, X2)P, and its 
breakup component, |^'bu(3^i5 a;2)P, respectively. The density along the diagonal xi = X2 
line now increases significantly, as in the full correlated calculations shown in Figs. 6(d) 
and 6(j). Evidently, the final state interaction is essential in the breakup processes. 



bu (a) correlatedr_ 




ct (fm) 

Figure 11. (Color online) The breakup probability as a function of time. The 
upper and the lower panels correspond to the correlated and the uncorrelated cases, 
respectively. The solid line shows the total breakup probability, while the dashed and 
the dot-dashed lines correspond to decompositions of the breakup probability into the 
(be) and (cc) processes, respectively. See the text for details. 
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The total breakup probability as well as the probability for each breakup process 
are shown in Fig. 10 as a function of time. They are defined as 

/oo 
c/xi(ix2 |*bu(a;l,X2,^)|^ (25) 
-oo 

= |«(^")(t)r + El/3i'^nt)P+ E \lt'\t)\^ (26) 



nj^O n<n' 

(n.n'^O) 



and 



Pbc(t) = r dx^dx, \^t''\^i,^2,t)\' = E (27) 

P,e(t)= rrfXiCiX2|vl>(^)(^l,^2,t)P= E l72^(^)l'> (28) 

J —oo ^ I 



n<n' 
(n,n'#0) 



which are denoted by the solid, the dashed, and the dot-dashed lines, respectively. The 
upper and the lower panels correspond to the correlated and the uncorrelated cases, 
respectively. Notice Pbb = for the uncorrelated case because of the orthogonalization 
of the wave functions. For the correlated case, although Pbb is finite, it is considerably 
small and is not shown in Fig. 10. By comparing the upper and the lower panels of 
Fig. 10, one finds that the total breakup probability is increased due to the pairing 
correlation. This is the case especially for Pec, that is, for the two- neutron emission 
process. 

One also sees that the (be) component first increases as a function of time while 
the increase of the (cc) component is somewhat delayed. For the correlated case, the 
(be) component eventually decreases and the (cc) component takes over. This is a 
manifestation of the dominance of a sequential mechanism in the two-neutron breakup, 
if one intends to use a terminology of perturbation theory. Nevertheless, the strong final 
state interaction makes the dineutron-like emission the main breakup process as shown 
in Fig. 6. 

In Ref. p2], we have pointed out that the properties of the two-body subsystem 
with a neutron and the core nucleus play a decisive role in the Coulomb breakup of ^^Li 
and ^He nuclei. This is because the external field is so weak for the Coulomb breakup 
that only one of the neutrons makes a transition to other (continuum) single-particle 
states. This corresponds to the (be) process in our example. The external field is much 
stronger for the nuclear breakup process, and the two-step process, or even higher step 
processes, play an important role. This can be seen in a large probability for the (cc) 
process. Pec, shown in the upper panel of Fig. 10. The effect of dineutron correlation 
can therefore be much easily seen in the nuclear breakup process as compared to the 
Coulomb breakup. A similar conclusion has been reached also in Ref. [51] (see also Ref. 



5. summary 



We have developed a simple schematic model for two-neutron halo nuclei, which still 
contains the essential features of physics of unstable nuclei. The model is based on a 
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three-body model in one spatial dimension. That is, the two valence neutrons move 
in a one dimensional mean field potential while interacting with each other via a two- 
body interaction. The two-body interaction scatters the two neutrons into many single- 
particle states with opposite parity, causing the strong dineutron correlation in the 
ground state. Applying this model to a weakly bound two-neutron halo nucleus, we 
have shown that the two-particle density for the ground state is indeed concentrated 
in the region of a;i ~ 0:2, as a consequence of the dineutron correlation. This model 
also allows detailed studies on several dynamical processes. We have solved the time- 
dependent two-particle Schrodinger equation in a non-perturbed way under the infiuence 
of a time-dependent external field, which can simulate the field generated by the reaction 
partner during a heavy-ion collision. We have shown that the main breakup process is 
an emission of dineutron, that is, the correlated neutron pair, in spite of the one-body 
nature of the external field and although the one neutron emission process is dominant at 
the early stage of time evolution. We have also shown that the pairing correlation, and 
thus the dineutron correlation, significantly enhances the breakup probability, especially 
for the two-neutron emission process. 

A schematic model such as that presented in this paper is useful to get deep insight 
into the physics behind. In addition to the nuclear breakup process studied in this 
paper, the one dimensional model for two-neutron halo nuclei can also be used e.g., in 
order to clarify the pair transfer reactions of exotic nuclei, whose dynamics has not yet 
been fully understood. A work in this direction is now in progress, and we will publish it 
in a separate paper [53] . Subbarrier fusion reaction of two-neutron halo nuclei is another 
interesting application of the one dimensional model. It would be straightforward to 
extend the one dimensional model for fusion of one-neutron nuclei [54] to two-neutron 
nuclei. Such study will shed light on the effect of irreversible process such as breakup 
on many-particle quantum tunneling. 
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Appendix A. Admixture of even- and odd-parity states 

In this appendix, we prove by using a two-level model that a and /3 in Eq. IHM 
have the same sign for the density-dependent pairing interaction. These coefficients 
are determined by diagonalizing a 2x2 matrix. 





with 



00 



A 




(A.2) 



—00 
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B 
C 

9{x) 



2e. 



9[l 



dx g{x) (f)eixY<Po{xy 
dxg{x) (j)o{x)^, 

1 



X _|_ ^(\x\-R)/a 

(see Eq. ([7])). The lower eigenvalue of this matrix is 
1 



A 



{A + C)-JiA-Cy + AB^ 



and the corresponding eigenvector is 



a 

/3 



2B 



(C- A) - J{C-Ay + AB^ 



(A.3) 
(A.4) 
(A.5) 

(A.6) 
(A.7) 



where N" is the normalization coefficient. Since both B and {C — A) — y {C — Ay + 
are negative, a and /3 have the same sign, as in the BCS ground state wave function in 
which all the configurations are superposed coherently. 
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